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Sponges are regarded as the most primitive 
multicellular animals, primarily due to their 
‘simple’ body plans and evidence from the fossil 
record, which together suggest they were the 
earliest lineage to diverge within the Metazoa. 
Sponges have attained a ‘multicellular’ grade of 
construction, with no development of tissues or 
organs, and their fossil record extends for at least 
600 million years. Due to their shared possession 
of unique choanocyte cells, the Porifera is gener- 
ally considered to be a monophyletic phylum. 
Recent ultrastructural and molecular evidence, 
however, suggests they may be a paraphyletic or 
even polyphyletic group of animals (Mackic & 
Singula, 1983; Reiswig & Mackie, 1983; 
Cavalier-Smith et al., 1996; Borchiellini et al., 
1998; Kruse et al., 1998). 


Phylum Porifera, as it is presently recognised, 
consists of three classes; Demospongiae, Hexact- 
inellida and Calcarea, differentiated primarily by 
differences in composition and geometry of 
skeletal components and cellular organisation of 
the soft parts. Hexactinellida and Demospongiae 
are characterised, in part, by the common 
possession of inorganic skeletons composed of 
siliceous spicules, although some demosponges 
(Dictyoceratida, Dendroceratida and 


Verongida), have skeletons composed of only 
proteinaceous (spongin) fibres and collagen fibrils. 
Demospongiae and Hexactinellida differ in that 
the former have spicules with one to four rays 
(monactine to tetractine), whereas the latter 
always have triactine or triaxial-derived 
(pentactinal and hexactinal) spicules. In contrast, 
Calcarea include sponges with calcium carbonate 
spicules in the form of calcite. Calcarea and Demo- 
spongiae both have representative species of 
‘coralline sponges’ possessing a calcareous 
aragonitic base, in addition to calcareous or 
siliceous spicules, respectively. These sponges, 
now referred to as ‘hypercalcified’, formerly 
comprised the Class Sclerospongiae, whereas it is 
now believed that this grade of construction has 
evolved independently in different lineages within 
the two classes (Vacelet, 1985; Reitner, 1992), 


Demospongiae and Calcarea have three major 
cellular layers. The first layer, the pinacoderm, 
lines all external surfaces of the sponge and is 
composed of a single layer of pinacocyte cells. 
The second layer is the choanoderm, composed 
of choanocytes which are the collared cells that 
draw water, and hence nutrition, into the sponge 
via the aquiferous canal system. Lastly, is the 
mesohyl, a proteinaceous matrix lying between 
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the pinacoderm and choanoderm, where the 
skeletal material is found with all other cell types. 


Morphologically, Hexactinellida are consider- 
ably different from Demospongiae and Calcarea, 
with syncytial cellular organisation. Instead of 
pinacocytes, these sponges have a syncytial 
surface dermal membrane which is contiguous 
with an inner trabecular membranc that drapes 
through the sponge interior. The mesolamella in 
hexactinellids is equivalent to the mesohyl in 
other sponges (Mackie & Singula, 1983). The 
mesolamella is composed of collagenous sheets 
which form a suspensory network for attachinent 
and support of trabecular tissues (Reiswig & 
Mackie, 1983). Hexactinellids do not possess a 
choanoderm as in the other two classes, but have 
a choanosyncytium composed of numerous 
collared bodies sharing a common nucleus and 
joined by stoloniferous cytoplasmic bridges. 
Hexactinellids also possess a unique secondary 
suspensory network that supports the collars of 
the collar bodies, which is not present in Demo- 
spongiae or Calcarea (Rciswig & Mackie, 1983). 
Because of these major differences Bergquist 
(1978) and Reiswig & Mackie (1983) proposed 
separate phylum and subphylum status, respect- 
ively, for the Hexactinellida. 


Two major hypotheses have been developed 
explaining the evolutionary relationships 
between the sponge classes (Fig. 1). The first 
suggests that Demospongiae and Hexactincllida 
are more closely related to each other than to 
Calcarea, based on their respective similarities in 
the chemical composition of spicules (Mohn, 
1984; Béger, 1988; Ax, 1996). Calcareous 
spicules are formed extracellularly by several 
sclerocytes (Ledger & Jones, 1977), whereas silic- 
eous spicules of Demospongiae originate intra- 
cellularly within a single sclerocyte, and those of 
Hexactinellida are formed intrasyncytially by a 
‘scleroblast mass’ containing many nuclei. The 
central axial filaments of spicules from Demo- 
spongiae and Hexactinellida differ in their cross 
sectional geometry (hexagonal and triangular vs. 
square, respectively), while calcareous spicules 
arc devoid of a central filament. These differ- 
ences in the axial filament could be indicative of 
differences in chemical composition (Reitner & 
Mehl, 1996). Most authors believe that spicules 
were derived independently in each of the three 
classes, and consequently are not useful as phylo- 
genetic character at the class level except to show 
that spiculogenesis is not homologous in each 
class. 
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FIG. 1. Diagrammatic representation of the two major 
hypotheses of relationship between the three 
poriferan classcs. A, Class Demospongiae and Class 
Hexactinellida are more closely related to each other 
than to the Calcarea based on chemical composition 
of the spicules (after Méhn, 1984; Béger, 1988; Ax, 
1996). B, Class Calcarea and Class Demospongiae 
(Subphylum Cellularia) are more closely related to 
each other than to the Hexactinellida (Subphylum 
Symplasma) due to their differing cellular condition 
(after Reiswig & Mackie, 1983), 


The second hypothesis proposes that Calcarea 
and Demospongiae are more closcly related to 
each other than to Hexactinellida, emphasising 
major differences in cellular condition between 
these groups. It was proposed that the Subphyla 
Symplasma (representing hexactinellids) and 
Cellularia (demosponges and calcareans) be 
erected to distinguish between these groups 
(Reiswig & Mackie, 1983). Reitner & Mehl 
(1996) proposed the term Pinacophora, stating 
that it was more appropriate than the term Cellularia 
when comparing sponges with other metazoans. 
They identified three apomorphies for the group 
which separated it from Hexactinellida: 1) the 
presence of a pinacoderm, 2) ball-shaped choano- 
cyte chambers in the adult, and 3) the ability to 
produce a calcareous (‘hypercalcified’) basal 
skeleton, as in the coralline sponges (Reitner & 
Mehl, 1996). 


Hexactinellida first appeared in the fossil record 
during the Late Proterozoic, approximately 
540my ago, whereas Calcarea and Demospongiae 
did not appear until some 50my later, in the Early 
Cambrian (Finks, 1970). It has since been 
proposed, however, that precursors to extant 
sponges were devoid of spicules and therefore 
would not have fossilised easily (Vacelet, 1985), 
which lends apparent support to the second hypo- 
thesis. Moreover, the known fossil history for 
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sponges has recently been challenged by Li et al. 
(1998), who claim to have found fossil demo- 
sponges 580my old. They suggest that demo- 
sponges were the first class to evolve, rather than 
hexactinellids (Li et al. 1998), supporting the first 
hypothesis. 


As demonstrated by Van Soest (1987), sponge 
classification at the higher levels is problematic 
due to a lack of synapomorphic characters and 
numerous assumed homoplasies. Characters such 
as spicule composition and cellular construction 
may be valid and truly indicative of phylogeny, 
but it ts difficult to determine which character 
should be given more weight, if at all. 


Recent molecular studies have been very 
helpful in providing additional characters to 
assist with phylogenetic reconstructions. In an 
attempt to elucidate the phylogenetic history for 
Porifera, various gene sequences have been 
explored including the small subunit of the 
ribosomal 18S gene (18S rRNA), heat shock 
protein 70 (Hsp70) and Protein Kinase C. 


Cavalier-Smith et al. (1996) analysed full- 
length 18S rRNA genes which yielded a 
paraphyletic Porifera. A Demospongiae- Hexact- 
inellida clade formed a sister group to a 
Calcarea-Ctenophora clade. This suggests that 
Demospongiae and Hexactinellida are more 
closely related to each other than they are to 
Calcarea. This hypothesis contradicts the proposal 
for subphylum status for Hexactinellida, as 
suggested by Reiswig & Mackie (1983) and 
Reitner & Mehl (1996). 


Analysis of the two protein coding genes, 
Hsp70 (Borchiellini et al., 1998) and Protein 
Kinase C (Kruse et al., 1998) produced a poly- 
phyletic and paraphyletic Porifera, respectively. 
Koziol et al. (1997) also examined the Hsp70 
gene, but considered it to be too conservative for 
resolution within the Porifera. Borchicllini et al. 
(1998), however, examined the Hsp70 gene using 
the first and second codon positions and found, 
with low bootstrap support, that Calcarea and 
Demospongiae formed a clade to the exclusion of 
Hexactinellida. Their results are controversial in 
that sponges were shown to be a group derived 
from other Metazoa, Cnidarians were hypoth- 
esised as being the first metazoans to diverge, 
followed by Ctenophora which formed a sister 
group to sponges. In some analyses, 
hexactinellids formed a clade with ctenophores, 
rather than with other sponge groups, while the 
18S rRNA gene showed ctenophores to be more 
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closely related to Calcarea (Cavalier-Smith et al., 
1996). Analyses of the Hsp70 gene always resulted 
in a Demospongiac-Calcarea clade, which supports 
the concept of the Subphylum Cellularia. 


Results from analysis of the Protein Kinase C 
gene Were more similar to the 18S rRNA data, 
with Calcarea forming a clade with the lower 
metazoans. These data also showed that hexact- 
inellids were the first to diverge from the meta- 
zoans, While demosponges formed a sister group 
to a calcarean-metazoan clade (Kruse et al, 
1998), Hexactinellida did not form a clade with 
demosponges, but instead formed a sister group 
to all other metazoans. Unfortunately, analyses 
for each gene yielded differing topologies, each 
with low bootstrap support, leaving the alleged 
phylogenetic relationships of sponges unresolved. 
Even though these studies yielded conflicting 
phylogenetic patterns, it is still possible that Por- 
ifera may not be monophyletic, as traditionally 
believed. 


To date, these are the only molecular studies 
which have included representatives from each of 
the three classes of sponges (although the only 
hexactinellid 18S rRNA sequence has not yet 
been made available in GenBank; West & 
Powers, 1993). The 18S rRNA gene is the most 
extensively studicd gene, with the largest 
database available for comparison. Universal 
primers have also been developed, making this 
gene relatively easy to obtain sequences in a short 
period of time. Although it has been suggested 
that the 18S rRNA gene docs not have enough 
signal to address the phylogenetic history of the 
lower metazoans, due to over saturation (Rodrigo 
ct al., 1994), it has been demonstrated that 
increased taxon sampling can assist in resolving 
phylogenetic relationships, principally by 
spreading homoplasic signal among a greater 
diversity of internal branches (Hillis, 1996). 


The aims of this study were to test the mono- 
phyly of Porifera by examining relationships 
between Demospongiae, Calcarea and Hexact- 
inellida. In addition to those six sequences 
already available in GenBank, we generated six 
additional full-length 18S rRNA sequences, 
strengthening inference from all available data. 
Phylogenetic trees inferred from Distance Matrix 
(DM), Maximum Likelihood (ML) and Maximum 
Parsimony (MP) methods were compared with 
hypotheses generated from other genes, utilising 
Constraint Analysis via the Kishino-Hascgawa 
test (Kishino & Hasegawa, 1989). 
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TABLE 1. List of sponges used in this study with 
corresponding museum voucher number and 
collection locality (identified by MK). 


‘i 3 Collecti 
Classification | Voucher number lo aa 
Class Demospongiae 

Subclass Tetractinomorpha p 

Vetulina stalactites : 

Schmidt, 1879 BMNH 1998.3.19.1| Caribbean Sea 
Acunthochaetetes wellsi ; Palau, 
Hartman & Goreau, 1975 |BMNH 1995.11.2.2 Micronesia 


; j 
Subclass Ceractinomorpha | 


ee 
acon ies MEBIN eee 
ia BMNH 1998.3.19.2| Red Sea 
Class Hexactinellida 

| Subclass Hexasterophora _| 
Syinpagella mex Schmidt, | HROM 003:00925 | Turks & Caicos 


Margaritella coelopty- 
chivides Schmidt, 1880 


Class Calcarea 


HBOM 003:00929 | Turks & Caicos 


Subclass Calcinea 


HBOM 


Leucetta sp. Bahamas 


27:X:96:3:305 


MATERIALS AND METHODS 


SAMPLE COLLECTION AND SELECTION. 
Specimens were collected by SCUBA and 
manned-submersible from various localities 
between 1989-1996 (Table 1). Immediately upon 
collection a_small piece of the sponge, approx- 
imately 3cm”, was removed from the interior with 
aclean scalpel blade to minimise surface epibiont 
contamination. This voucher was either frozen at 
-20°C or diced as finely with a sterile razor blade 
and immediatcly placed in Guanidium Chloride 
(GnCl) Buffer [6M GnCl, 5% Tween 20, 0.5% 
Triton X-100 in 1L Tris-EDTA buffer, pH 8.0 
(100mM Tris, 30mM EDTA)] for stable storage 
of the lysed cells and DNA. A representative 
piecc of the sponge was placed in 70% ethanol for 
subsequent taxonomic identification. Voucher 
specimens of all taxa were deposited at The 
Natural History Museum, London (BMNH), the 
Harbor Branch Oceanographic Institution 
Muscum, Fort Pierce, Florida(HBOM) and in the 
personal collection (MKB) of MK. 


EXPERIMENTAL PROCEDURES. Sponge cell 
buftered lysate was diluted four-fold with auto- 
claved analar H,O, and frozen specimens were 
thawed and ground with a micropestle to form a 


slurry prior to DNA extraction, using standard 
phenol:chloroform extraction procedures 
(Sambrook et al., 1989). Full-length (approx- 
imately 1300bps) 18S rDNA was amplified with 
the forward and reverse primers 18Sf20 and 
18Sr21, respectively (McInerney et al., in press). 
PCR conditions for 1.6ng/ul DNA in 50ul 
reactions were: initial denaturation at 94°C for 
5mins, followed by 35 cycles of denaturation at 
94°C for Imin, annealing at 55°C for 1min, and 
extension at 72°C for Imin. The product was 
electrophoresed on an 0.8% agarose gel stained 
with lug/ul ethidium bromide to check band 
size, and then purified from the gel with the 
Qiaex 11 PCR purification kit (Qiagen Ltd, UK), 
following manufacturers instructions. In addition 
to the two PCR primers, eight internal primers 
were used to sequence both chains automatically, 
utilising the dideoxy chain termination method 
(Sanger et al., 1977) (forward primers 377F 
CCGGAGARGGAGCCTGA, 577F GCCAGC 
MGCCGCGGT, 1262F GGTGGTCGATG 
GCCG and 1510F CAGGT CTGTGATGCC and 
their complementary reverse primers called 
377R, 577R, 1262R and 1510R). Each 
contiguous sequence fragment was replicated 
with at least one overlapping fragment (Amersham 
Cycle Sequencing Kit), Ambiguous nucleotide 
positions were coded according to the Inter- 
national Union of Pure and Applied Chemistry 
(IUPAC) nomenclature. Sequences were managed 
utilising Sequencher 3.0 software (Gene Codes 
Corporation, 1995). New sequences were 
deposited in the GenBank sequence repository 
(http://www?.ncbi.nlm.nih.gov) under accession 
numbers Verulina stalactites (AF084236), 
Acanthochaetetes wellsi (AF084237), Clathria 
(Thalysias) reinwardti (AF084238), Negombata 
corticata (AF084239), Sympagella nux 
(AJ224123), Margaritella coeloptychioides 
(AJ224124) and Leucetta sp. (AF084240). 


PHYLOGENETIC RECONSTRUCTION. To 
test the monophyly of Porifera, sequences from 
the lowcr metazoans were included in the 
analysis: Phylum Porifera - Tetilla japonica 
(D15067), Microciona prolifera (L10825) 
(which Hooper (1996) referred to the subgenus 
Clathria (Clathria) based on taxonomic 
re-evaluation of type material, and we assume 
that GenBank L10825 belongs to this specics), 
Axinella polypoides (U43190), Clathrina 
cerebrum (U42452), Seypha ciliata (L10827) 
(which belongs to Sycon, follwing Dendy & 
Row, 1913; Gert Woerheide, pers.comm.), Sycon 
calcaravis (D15066), Phylum Placozoa - 
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Trichoplax adhaerens (L10828), Phylum 
Ctenophora - Beroe cucumis (D15068), 
Mnemiopsis leidyi (L10826); Phylum Cnidaria - 
Anemonia sulcata (X53498), Tripedalia 
cystophora (L10829). Representatives of each of 
the major fungal groups were chosen as the out- 
group taxon, due to their inferred position as a 
sister group to Mctazoa, according to rRNA data 
(Wainright et al., 1993) and protein data (Baldauf 
& Palmer, 1993): Fungi - Aureobasidium 
pullulans (M55639), Saccharomyces cerevisiae 
(Z75578), Athelia bombacina (M55638), 
Blastocladiella emersonii (X54264). Taxa selcc- 
ted for this study were retrieved from a sccondary 
structure alignment maintained on the Ribosomal 
Database Project (RDP) database (http://rdpwww. 
life.uiuc.edu/ index2.html; Maidak et al., 1996). 
The profile alignment option of ClustalW was 
then used to combine the two alignments (Higgins 
& Sharp, 1988). Sequcnces were aligned using 
ClustalW 1.7 and then modified by eye in the 
Genetic Data Environment, GDE, (Smith et al., 
1994) on a SUN workstation. A conservative 
approach was used for alignment. Only those 
positions whose alignment was ambiguous were 
choscn for analysis, thus ruling out a significant 
number of potential positions. Approximately 
1300bps were sequeneed. The new alignment 
(including other sponges, inetazoans and fungi), 
was 2590bps in length, and after removal of 
ambigous positions the resulting alignment was 
987bps long. 


Phylogenetic hypothescs werc constructed, and 
sequence statistics werc cvaluated, using PAUP* 
4.0.0d64 test version (Swofford, in press). The 
likclihood ratio test statistic was used to evaluate 
which evolution model fit the data best 
(Goldman, 1993). This was achicved by first 
constructing a neighbour joining tree calculated 
by the Jukes & Cantor (1969) method. Using this 
guide tree and maximum likelihood criteria. the 
transition-transversion ratio, base composition 
and proportion of invariable sites were estimated 
using the Newton-Raphson method implemented 
in Paup*4.0. Each of these variables were 
calculated separatcly, and then entercd into the 
model to calculate the next variable. This proccss 
was repeated until the best maximum likelihood 
value was rcached. The chosen evolutionary 
model was that which yielded the optimal likeli- 
hood value without compromising model 
complexity. The reliability of internal branches 
was evaluated by the bootstrap resampling 
method (Felsenstein, 1985). In each analysis, 100 
iterations were carried out for each optimality 


criterion. A 50% majority rule consensus tree 
was inferred from the resulting bootstrap 
partition table. 


The inferred phylogenetic relationships deriv- 
cd trom the new I8S rRNA sequences were 
compared with the major phylogenetic hypo- 
theses derived from morphological characters 
and other gencs. Using MacClade (Maddison & 
Maddison, 1992), we constructed trees representing 
competing hypotheses. The constrained trees were 
examined using the Kishino-Hasegawa test 
(Kishino & Iascgawa, 1989), lor both ML and 
MP methods. 


ABBREVIATIONS. DM, Distance Matrix; ML, 
Maximum Likelihood: MP, Maximum 
Parsimony methods. 


RESULTS 


Prior to bootstrapping, MP, DM and ML 
methods yielded trees with virtually identical 
topologies, except for the branch arrangement 
within Demospongiae which was largely unres- 
olved. For the DM method, the F84 model was 
used (Felsenstein, 1984), with the minimum 
evolution objective function. The proportion of 
sites assumed to be invariable = 0.686332. For 
the MP method, (parsimony informative sites = 
155), the characters were treated as unordered 
and equally weighted. The ML analysis was 
conducted with a Two-type substitution model 
with an estimated transition/transversion ratio of 
1.553481, estimated base frequencies of A = 
0.270333, C = 0.206219, G = 0.258425 and T = 
0.265023, and an estimated proportion of 
invariable sites = 0.686760 (with equal rates of 
variation for all sites). Using the MP method as 
representative, two equally parsimonious trees 
with minimal differcnces in tree topology were 
recovered (C1 = 0.669, RI = 0.723, tree length = 
465, total number of characters = 987, parsimony 
informative characters = 155) (Fig. 2). The only 


- difference in branch arrangement between the 


two MP trecs was the position of Acantho- 
chaetetes within the Demospongiae. In Tree A 
(Fig. 2), Acanthochaetetes is the earliest 
separation within Tetractinomorpha, whereas in 
Tree B (Fig. 2), Acanthochaetetes is the carliest 
separation of the Ccractinomorpha clade. The 
three sponge classes are monophyletic in both 
trecs, with Demospongiae as sister taxon to 
Hexactinellida, while Calcarea formed a sister 
taxon to Ctenophora. 

Within Demospongiae, Subclass Tetractino- 
morpha is represented by Verulina and Tetilla, 


Acanthochaetetes 
Vetulina 
Tetilla 
Axinella 
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| Tetractinomorpha 


Clathria (Thalysias) : 
Clathria (Clathria) | Ceractinomorpha Porifera 
Negombata 
Sympagella | Hexasterophora 
Margaritella i 
Sycon calcaravis | Calcaronea 
Sycon ciliata A 
à Clathrina | Calcinea 
Leucetta 
Mnemiopsis | Ctenophora 
Beroe 
Anemonia | Cnidaria 
Tripedalia 
Trichoplax 1 Placozoa 
Aureobasidium 
Saccharomyces Fungi 
Athelia 
Blastocladeilla 
‘Acanthochaetetes f 
Clathria (Thalysias) Ceractinomorpha 
Clathria (Clathria) 
Negombata 
Vetulina Tetractinomorpha 
Tetilla Porifera 
Axinella home 
Sympagella | Hexasterop 
Margaritella EN 
Icaravis 
B Sycon ciliata : 
Clathrina | Calcinea 
Leucetta 
Mnemiopsis | Ctenophora 
Baos i 
TE | Cnidaria 
Tripedalia 
į Placozoa 
Aureobasidium 
Saccharomyces Fungi 


FIG. 2. Two most parsimonious trees derived from molecular data analysed using the assumptions of MP (CI 


lastocladiella 


0.669 , RI=0.723, tree length = 465, total number of characters = 987, parsimony informative characters = 155). 
A, Acanthochaetetes is the earliest taxon to diverge from the Tetractinomorpha. B, Acanthochaetetes is the 


earliest taxon to diverge from the Ceractinomorpha. 


which form a common clade. Axinella was also 
located within the tetractinomorphan clade. The 
family Axinellidae was formerly considered to be 
a tetractinomorph on the basis of their repro- 
ductive strategies, but Van Soest et al. (1990) 
suggested that, based on morphology, there was 
more support for their placement among 
Ceractinomorpha. Our dataset, however, places 
Axinellidae (represented by Axinella), in Tetrac- 
tinomorpha. Acanthochaetetes is also 
traditionally recognised as a tetractinomorph 


sponge, but this position is unstable, grouping 
with tetractinomorphan Tree A, and with 
Ceractinomorpha in Tree B. Ceractinomorpha 
were represented by two species of Clathria, and 
Negombata, where the latter taxon has had a 
questionable affinity with switches between 
Tetractinomorpha (Bergquist, 1978) and 
Ceractinomorpha (Topsent, 1922). Kelly-Borges 
& Vacelet (1995) suggested that based on morph- 
ology, chemistry and reproduction, Negombata 
has a closer affinity to Ceractinomorpha, Order 
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TABLE 2. Results of the K-H test for alternatives to the MP/ML tree. 
Tree number | is the best tree in both analyses. Tree 2 is a 
constrained tree that retains a clade containing all of the sponge taxa. 
Tree 3 is a constrained tree that places Calcarca with Demosponges 
to the cxclusion of all other taxa. Tree 4 places Calcarea with 
Hexactinellida. The second column gives the ML score (negative 
log-likelihood) for cach tree. The third column gives the difference 
in Log-likelihood between alternatives. The fourth column 
indicated whether or not an alternative is significantly worse than 
the optimal tree (P<0.05). The fifth column gives the tree length of 
the various trees. Column number six gives the difference in tree 
length between alternative trees. The last column indicates the 
significance of the result (an asterisk indicates a P-value < 0.05). 


trees relating most Demospongiae 
taxa received less than 50% 
bootstrap support. The Negombata - 
Clathria clade was the only 
relationship that remained intact and 
formed a monophyletic group, with 
bootstrap values of 98 for DM, 100 
for MP and 99 for ML, yielding 
further support to the hypothesis that 
Negombata is more closely related to 
Poecilosclerida than to the 
tetractinomorph Hadromerida. 


E E ed T | Support for a Demospongiae- 
' Tree = aa] Hexactinellida clade was relatively 
-Ln L Di -In L | pr Length | Diff. pe | high (DM 89%, MP 81% and ML 

1 | 3960.28210 | (besi) {| 466 | (best) | _ 70%), whilc Calcarea formed a 

2 | 3964.75361 | 447151 | 0.4971 | 469 3 0.3176 | monophyletic group with bootstrap 

3 | 4128.96432 | 168.68222 | <0.0001* 514 i 48 |<0.0001*| support of 100 for the DM and MP 

4 | 412238473 | 162.10263 | <0.0001* sio | 44 |[<o.0001%} trees and 98 for the ML tree. 


Poecilosclerida. Our data lend support to this 
hypothesis. 


There were only two representative taxa for 
Hexactinellida, each belonging to the same sub- 
class, Ilexasterophora, forming a monophyletic 
clade in all three analyses. Hexactinellida consist- 
ently form a sister group with Demospongiae. 


The topology of trces obtained from the three 
methods of analysis varied only slightly from that 
of the DM tree, recovering a topology which 
suggests that the calcareous Subelass Calearonea 
arose from within the Subclass Calcinea (not 
shown). Both MP and ML methods yielded a 
topology with the two calcareous subclasses posi- 
tioned as sister taxa. Our data show the 
Subclasses Calcaronea and Calcinea to be valid 
groupings, but at Icast three taxa per group are 
needed to infer relationships. Calcarea form a 
sister group relationship with Ctenophora. This 
clade is hypothesised to be derived from the 
diploblastic animals (Cnidaria and Placozoa). 


A total of 100 bootstrap replicates were constr- 
ucted for each of the three methods of analysis; 
DM, MP and ML (Fig. 3). Each of the three 
classes is consistently recovered, with high 
bootstrap support as a monophyletic elass, Hexact- 
inellida formed a clade with 100% bootstrap 
support for each method, and Demospongiac 
formed a monophyletic group with bootstrap 
values of 100, 99 and 88 for DM, MP and ML 
methods, respectively. There were no 
representatives of the Subclass Homo- 
scleromorpha in these analyses due to difficultics 
we encountered during PCR. The topology of 


Calcarea consistently formed a sister 

taxon with Ctenophora, although 
there was only low bootstrap support of 64, 55 
and 56, for DM, MP and ML methods, 
respectively. 


In order to test alternative hypotheses on 
possible relationships between the sponge classes, 
the Kishino-Hasegawa (K-H) test was employed 
for both the MP and ML methods. Constraint 
trees were designed and compared against the 
optimal tree inferred from our data, which 
showed Calearea was more closely related to 
Ctenophora, and Demospongiae and Hexactin- 
ellida formed a sister-group. The optimal tree for 
each constraint analysis was found using the MP 
and ML methods, and then the resulting tree was 
compared with the optimal, unconstrained tree to 
test for significant differences (Table 2). 


For the K-H test using ML and MP methods, 
three constrained trees were designed and 
tested against the optimal tree derived from the 
molecular data (Tree 1): a monophyletic clade 
containing all three classes (Tree 2); a clade 
containing the I lexactinellida as a sister taxon to a 
Demospongiae-Calcarea clade (Tree 3); and a 
clade containing the Demospongiae as a sister 
taxon to a Hexactinellida-Calcarea clade (Tree 
4). All alternative trees were rejected when 
compared to the aptimal tree, with the excep- 
tion of a monophyletic Porifera (Tree 2). 


The optimal tree from ML analysis has a log- 
likelihood value of -3960.28210, while the best 
tree which retained a monophyletic sponge clade 
had a log-likelihood score of -3964.75361. A 
two-tailed T-test did not reject this hypothesis as 
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FIG. 3. Phylogenetic consensus tree for DM, MP and ML with bootstrap values corresponding to each method. 
50% majority rule consensus trees are displayed. DM values are presented in normal text, MP values are in bold 


and ML bootstrap values are in italics. 


being significantly worse than the ML tree. 
Similarly, using the MP analysis method, the most 
parsimonious tree was only three steps shorter 
than a tree which contained a monophyletic 
sponge clade, and the Log-likelihood of a mono- 
phyletic tree was only 0.1% less likely than the 
optimal tree. In both instances, the alternative 
hypothesis was not significantly worse than the 
optimal tree (P<0.05). Given that we cannot reject 
the alternative hypothesis, it would be unwise to 
suggest that the phylum is not monophyletic. 


DISCUSSION 


Bootstrap results indicate a polyphyletic 
arrangement for Phylum Porifera, and support the 
theory that siliceous sponges evolved separately 
from calcareous sponges. This arrangement has 
previously been suggested in the literature, 
supported by data from 18S rRNA (Kobayashi et 
al., 1993; Cavalier-Smith et al., 1996), 28S rRNA 


(Lafay et al., 1992) and Protein Kinase C genes 
(Kruse et al., 1996). In these earlier studies, the 
sponge classes and subclasses were not exten- 
sively represented in analyses, whereas our data 
doubles the number of sponge 18S rRNA 
sequences analysed, and yet yields the same 
topology. Our results also show that calcareous 
sponges are derived from other metazoans, which 
are generally considered to have evolved later 
than the siliceous sponges, based on morpho- 
logical data and the fossil record. 


The poriferan classes have several apparent 
apomorphies, but upon closer examination these 
may be convergent characters. Differences in 
spicule geometry and spiculogenesis are obvious 
between calcareous and siliceous sponges, but 
possibly also between the Hexactinellida and 
Demospongiae. There are notable differences in 
choanocyte size, shape and arrangement within 
choanocyte chambers at the class and order 
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levels. While species of Demospongiae and 
Calcarea share similar cellular construction and 
mesohyl characteristics, when compared to 
Hexactinellida, Calcarca lack the collagen and 
proteinaceous fibre development found in Demo- 
spongiae. Larval morphology is also very 
different within and between the three classes. It 
is possible that convergent evolution has taken 
place between these three classes due to a 
common benthic, filter-feeding lifestyle. 


lt is also quite difficult to explain why Calcarea 
are consistently grouped with the Ctenophora 
according to molecular data, even with low 
bootstrap support. Morphologically, Ctenophora 
are much more complex than Porifera, with true 
lissue structure such as mesodermal muscle, 
gonoducts, and an anal pore; however increased 
complexity does not necessarily equate with a 
more derived evolutionary state. Itis conceivable 
that Calcarea have secondarily lost characters 
that are present in ctenophores. Cavalier-Smith et 
al. (1996) suggest that calcarcous spicules found 
in Calcarea may be homologous to those found in 
Cnidaria. A loss of spicules could have occurred 
in Ctenophora. They also suggest that larvae found 
in the subclass Calcinea (Calcarea) arc morph- 
ologically much more similar to those of other 
animals than they are to other sponges. 
Additionally, the position of ctenophores with 
respect to other diploblasts, based on molecular 
data, rests on analysis of only two representatives 
of the phylum. 


Although we failed to confirm whether the 
Porifera was a monophyletic taxon using our 
expanded molecular dataset, the hypothesis of 
poriferan monophyly is not rejected based on 
molecular data. The exact placement of 
calcareous sponges is problematic and requires 
further empirical support from sequenccs derived 
from a greater diversity of sponges and cteno- 
phorans. Surprisingly, when ctenophores were 
removed from analyses (data not shown), 
Calcarea remained a sister taxon to other 
diploblasts, and not with siliceous sponges 
(Demospongiae, Hexactinellida). Although these 
results are congruent with previous data, they are 
difficult to explain on the basis of ‘classical’ 
morphological characters and the paleontol- 
ogical record. 


Using the bootstrap resampling method, we 
have shown that a partition separating the Classes 
Hexactinellida and Demospongiae from all other 
taxa is reasonably well supported. Although 
bootstraps are not statistically high, any 


alternative arrangement receives very little 
support. This arrangement does, however, refute 
the hypothesis that Hexactinellida merit phylum 
or subphylum status. The conclusions of this 
study do not rule out the monophyletic nature of 
the phylum. Until data are gathered that can yield 
high confidence levels for a polyphyletic 
phylum. the possibility of sponges being mono- 
phyletic cannot be dismissed. 
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APPROACH TO THE PHYLOGENY OF 
AXINELLIDAE (PORIFERA: DEMO- 
SPONGIAE) USING MORPHOLOGICAL AND 
MOLECULAR DATA. Afemoirs of the Queensland 
Museum 44: 43, 1999:- A set of 27 specics of marine 
sponges of the Axinellidae and related families was 
selected with the aim of testing the monophyly of 
Axinellidae and investigating their phylogenctic 
relationships using cladistic methods. Partial 28S 
rDNA sequenccs, including the D3 domain, and 
traditional morphological characters were used 
independently to construct phylogenetic trees. 
Alignment of the sequences using the appropriate 
model of secondary structure of the RNA was 
compared to that produced by the ClustalW. The 
alignment using secondary structure constraints 
produced a better estimate of the phylogeny and was 
demonstrated to be an effective and objective method. 
The results from the analyses of the molecular and 
morphological data sets are not fully congruent: the 
morphological data suggest that Axinellidae is 
monophyletic, however the molecular data suggest 


that it is not monophyletic. In both cases the sampled 
members of the family arc closely related to those of 
Halichondriidae and Dictyonellidae. Tests of 
heterogeneity (reciprocal T-PTP and partition 
homogeneity test) shown that the data partitions are 
heterogeneous, which could be due to sampling errors 
(in either data set) or differences in the underlying 
phylogenies, and therefore data were not combined in 
a single analysis. O Porifera, Axinellidae, secondary 
structure, D3 domain 28S ribosomal DNA, phylogeny. 
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